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Abstract. Spontaneous creation of electron-positron pairs out of the vacuum due 
to a strong electric field is a spectacular manifestation of the relativistic energy- 
momentum relation for the Dirac fermions. This fundamental prediction of Quantum 
Electrodynamics (QED) has not yet been confirmed experimentally as the generation of 
a sufficiently strong electric field extending over a large enough space-time volume still 
presents a challenge. Surprisingly, distant areas of physics may help us to circumvent 
this difficulty. In condensed matter and solid state physics (areas commonly considered 
as low energy physics), one usually deals with quasi-particlcs instead of real electrons 
and positrons. Since their mass gap can often be freely tuned, it is much easier to 
create these light quasi-particles by an analogue of the Sauter-Schwinger effect. This 
motivates our proposal of a quantum simulator in which excitations of ultra-cold atoms 
moving in a bichromatic optical lattice represent particles and antiparticles (holes) 
satisfying a discretized version of the Dirac equation together with fermionic anti- 
commutation relations. Using the language of second quantization, we are able to 
construct an analogue of the spontaneous pair creation which can be realized in an 
(almost) table-top experiment. 



Optical lattice quantum simulator for QED in strong external fields 



2 



1. Introduction 

Spontaneous creation of electron-positron (fermion-antifermion) pairs from vacuum 
under specific external conditions is a direct manifestation of the relativistic energy- 
momentum relation for the Dirac particles. The most prominent realization of this effect 
is creation and separation of an electron (e~) and a positron (e^) in presence of a strong 
electric field, derivable as a phenomenon within Quantum Electrodynamics (QED). The 
adiabatic character of the spontaneous pair creation allows for the interpretation in 
which the particle is slowly pulled from the otherwise unobservable Dirac sea while 
the hole in the sea appears as an antiparticle. Both come into being on the cost 
of the external fields. Unfortunately, electrons and positrons, the lightest fermions 
satisfying the Dirac equation [1], still defend themselves from being exposed in that 
way. Generation of a strong enough electric field which is able to deliver the minimal 
energy 2m^c^ = 1.022 MeV needed to create a pair from vacuum still appears to be an 
experimental challenge. A natural source of a strong and localized electric field, the 
atomic nucleus, would need to carry a charge of at least +173 e (slightly depending on 
its predicted size, see e.g. [2]) what is about 50 unit charges above the heaviest (and 
unstable) nuclei which have ever been observed in the laboratory. In the early 1980s, 
there have been serious experimental attempts [3] to collide beams of fully ionized 
uranium atoms f/^^+ or similar ions in order to create a sufficiently long-lived charge 
concentration of around +184 e but they were not successful in this regard. Recent 
developments - e.g., in the field of strong lasers [4] or the current extension of GSI in 
Darmstadt - have again renewed interest in spontaneous pair creation. However, we 
are clearly not yet in the position of creating electric fields of sufficient strength. 

Quite surprisingly, help may come from distant areas of physics: condensed matter 
and solid state physics - areas commonly considered as low energy physics. Since the 
energy scale is determined by the mass of the particles under consideration, the electron 
mass can set a too high barrier for electrons and positrons while the analogous gap can 
be much lower for light quasi-particles whose masses can be tuned in experiments. 
This motivates our proposal of a quantum simulator in which excitations of ultra- 
cold atoms moving in a regular optical lattice will represent particles and antiparticles 
(holes) satisfying a discretized version of the Dirac equation together with fermionic anti- 
commutation relations. Applying the language of second quantization, we construct an 
analogue of the spontaneous pair creation which can be realized in an (almost) table-top 
experiment. 

To additionally motivate the need of a quantum simulator, we mention some open 
problems still present in theory and experiment related to supercritical fields of QED. 
The simplest setting in which the spontaneous pair creation should occur is the case of a 
constant electric field E, well known in the literature as the Schwinger effect or Sauter- 
Schwinger effect [5, 6, 7]. For nonzero values of the electric field -E > one should 
observe spontaneously generated pairs of particles and antiparticles with probability 
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(per unit time and volume) given by 

p...- ~ exp |- ^ ^ j = exp |- ^1 , (1) 

where £"3 = M'^c^ / (Hq) is the critical field strength determined by the elementary charge 
q and the mass M of an electron (or positron). Besides the aforementioned experimental 
difficulties, the above expression for Pg+e- is non-perturbative in qE and does not permit 
any expansion in the field strength E nor in the coupling constant (or charge) q, e.g. via 
a finite set of Feynman diagrams. Thus, apart from the constant field case, only very 
simple field configurations, where the electric field either depends on time E{t) or on one 
spatial coordinate such as E{x), have been treated analytically so far [8]. Consequently, 
our theoretical understanding of various aspects of this effect under general conditions 
is still quite limited. For example, recently it has been found that the occurrence of 
two different frequency scales in a time- dependent field E{t) can induce drastic changes 
in the (momentum dependent) pair creation probability [9, 10]. Moreover, the impact 
of interactions between the electron and the positron of the created pair, as well as 
between them and other electrons/positrons or photons is still not fully understood. 
This ignorance is unsatisfactory not only from a theory point of view but also in view 
of planned experiments with field strengths not too far below the critical field strength 
Es and thus capable of probing this effect experimentally [4]. 

The proposed quantum simulator will reproduce the quantum many-particle 
Hamiltonian describing electrons and positrons in strong electric fields and should 
thereby reproduce the Sauter-Schwinger effect. This will facilitate investigation of space- 
time dependent electric fields such as E{t, x) and also provide new insight into the role 
of interactions which may be incorporated into the simulator. 

It should be stressed here that our proposal goes beyond the simulation of 
the (classical or first-quantized) Dirac equation on the single-particle level, see, e.g., 
[11, 12, 13, 14, 15, 16, 17], but aims at the full quantum many-particle Hamiltonian. 
A correct description of many-body effects such as particle-hole creation (including 
the impact of interactions) requires creation and annihilation operators in second 
quantization. There are some proposals for the second-quantized Dirac Hamiltonian 
[18, 19, 20, 21, 22, 23] but they consider scenarios which are more involved than the 
set-up discussed here and aim at different models and effects. Similarly, the recent 
observation of Klein tunneling in graphene [24] deals with massless Dirac particles - but 
the mass gap is crucial for the non-perturbative Sauter-Schwinger effect, cf. Eq. (1). 
Furthermore, graphene offers far less flexibility than optical lattices regarding the 
experimental options for changing the relevant parameters or single-site and single- 
particle addressability, etc. 

2. Spontaneous pair creation in supercritical external fields 

We consider the Dirac equation [1] describing electrons/positrons propagating in an 
electromagnetic vector potential which are described by the spinor wave-function ^ 
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{h = c=i) 

r{id^-qA^)m-Mm = {}. (2) 

For simplicity, we consider 1+1 dimensions (/i = 0, 1) where the Dirac matrices 
satisfying the Chfford algebra {7'^,7'^} = It]^^ can be represented by Pauli matrices 
7*^ = (Ta and 7^ = —iai. Furthermore, we can choose the gauge qA^ = $ and Ai = 
(because in one spatial dimension, there is no magnetic field). In one spatial dimension, 
there is also no spin, hence the wave function has only two components = (\E'^, \E'^). 
As a result, the Dirac equation simplifies to 

idt^^it, x)=n ^(t, x) = {-ia2d^ + Mcxg + x). (3) 

If is negative and vanishes at infinity sufficiently fast the spectrum of the Dirac 
Hamiltonian "H consists of two continua (—00, — M] U [M, 00) and a discrete set of bound 
states En lying in the gap (— M, M). 




Figure 1. Typical spectrum of a Dirac Hamiltonian (left) and its dependence on the 
strength of an attractive potential $ (right). At the critical value $ = <I>cr the lowest 
bound-state Eq (red solid line) turns to a resonance in the negative continuum (red 
dotted line). Next bound states iJi, £'2, ... follow for larger values of $. 

The bound-state energies En depend continuously on the parameters of the 
potential. In particular, as the strength of the negative potential |$| increases each 
En \ —M. Already at a finite value called critical, the lowest lying bound 
state Eq reaches the negative continuous spectrum associated with the interpretation of 
antiparticles, i.e. £"0! = —M. For supercritical strength of the potential |$| > |$cr| 
the bound state - corresponding to a real pole in the resolvent of % (or in the scattering 
operator) - turns to a resonance (complex pole) with Re(i?o) < — M (see Fig. 1). 

Imagine now a time-dependent process in which $(t) is slowly varied between the 
sub- and supercritical regimes as in Fig. 2. In agreement with the adiabatic theorem, 
in the subcritical phase, the quantum state of the system follows the eigenstate in 
which it is initially prepared. As the supercritical phase begins the gap closes and 
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continuum 




Figure 2. Spectrum of the Dirac Hamiltonian T-L{t) in presence of a slowly varying 
potential In the middle, a supercritical phase: the lowest bound-state (red solid 

line) enters the negative continuum and turns to a resonance (red dotted line). 

the adiabatic theorem breaks down [25, 26]. The system follows then a resonance 
which is spectrally represented by a wave packet with position and width varying in 
time. Such wave packets inevitably decay in the lower continuum, trapping a big 
part of the wave function. Therefore, during the switch-off phase, when the potential 
$(t) becomes subcritical again, only a small part of the wave function follows the re- 
appearing eigenstate. Mathematically, there exists a non-vanishing matrix element of 
the scattering operator between the positive (-I-) and negative (— ) continuum, ^ 0, 
which tends to one in the "adiabatic limit". In order to avoid interpret ational problems 
we need to leave the one-particle picture at this stage and switch to the many-particle 
description. In the language of second quantization, the discussed process is described 
by the scattering operator S acting in a Fock space and it can be determined from the 
one-particle counterpart [27]. It involves dynamical and spontaneous creation of pairs 
as well as annihilation and scattering of already present particles and antiparticles. In 
the "adiabatic limit", the dynamical pair production goes to zero such that only the 
spontaneous process remains. Therefore, for processes (as described above) starting 
from a vacuum state and running through a supercritical phase we obtain [25, 28] 

— ^ S\n) = a^b^Q) (4) 

while for subcritical processes \Q) — ?■ in agreement with the adiabatic theorem. 
This phenomenon is called spontaneous pair creation, as opposed to the dynamical pair 
creation, since it is related to the spontaneous decay of a time- dependent ground state 
|fi(t)) (in the so called Furry picture) during the supercritical phase [29, 28]. 

3. Discretized quantum Dirac field 

In this section we make the first step towards the quantum simulator of a quantum Dirac 
field in optical lattices and discretize the theory by introducing a regular lattice in space. 
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We will argue that the phenomena of strong external fields like the supercriticality and 
spontaneous pair creation discussed above will survive this operation. 
The Hamiltonian for the classical Dirac field reads 

H = j dx ¥{-i(j2d^ + Mas + • (5) 

We introduce a regular grid (lattice) Xn = n ■ i with a positive grid (lattice) constant 
i and integers n G Z. The discretization of the wave function \&n(^) '■= V^^it, Xn), 
defined now at the grid points a;„, gives rise to a discretized derivative y/i dx^{t, Xn) — )■ 
[^n+i{t) — \E'.„_i(t)]/(2£) and to a discretized potential := $(x„). Finally, replacing 
the x-integral by a sum, we obtain 



-^i^n+l - ^n-l) + Maa^n + $n*r 



(6) 



In order to obtain the quantum many-body Hamiltonian, we quantize the discretized 
Dirac field operators via the fermionic anti-commutation relations 

{Kit), [K{t)n = SnmS'-^ , {nit),^im = . (7) 

Substituting = a„ and = 6.„ the discretized many-particle Hamiltonian obtains 
the form 

H=YiYl [^l+i^n - hW+i + h.c] + [($„ + M)didn + ($n " M)6],6„ . (8) 

n n 

The first term describes jumping between the neighboring grid points while the 
remaining two terms can be treated as a combination of external potentials. Due to 
the specific form of the jumping, the lattice splits into two disconnected sub-lattices: 
(A) containing and b2n+i and (B) containing a2n+i and b2n with integers n. Since 
the two sub-lattices behave basically in the same way, it is sufficient to consider only 
one of them, say A. With a re-definition of the local phases via a2n {—^)^CL2n and 
^2n+i (— l)"'''^^2n+i5 "wc obtain (for sub- lattice A) 

^ = ~ ^ X [^2n+l«2n + ^L-l«2n + h.C. + (9) 

n 

n 

Identifying C2n = fl2n and C2„+i = &2n.+i; this takes the form of the well known Fermi- 
Hubbard Hamiltonian for a one-dimensional lattice 

^ = X [^n+lCn + 4cn+l] + ^ ^"Cj,C„ , (10) 

n n 

with hopping rate J = l/£ and on-site potentials = + (— 1)"M. This Hamiltonian 
will be the starting point for the design of the optical lattice analogy. 

Alternatively, using a more abstract language of modern quantum field theory, the 
Hamiltonian H, being an operator acting on the Fock space J-", can be directly obtained 
by implementation of the discretized single-particle Hamiltonian 

^d^n = -^a2(^n+l - ^n-l) + McTg^n + $„^„ (11) 
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acting in the discretized Hilbert space S^d = (-^^(^))^ (which is the discretization of 
= L^(M)^) as a self-adjoint operator in the Fock space J-" according to 

H = Y, ^* Wn)^(/n) = ^ 5^ [Kll^l - ^V^Ul + h.C." 

(12) 

+ ^[($„ + M)^^t^i + (<|.„-M)^2t^2 

n 

where the second quantized discretized Dirac field operator 

^^=^"(/„) = j dx^^{x)f:{x) (13) 

satisfies the above anti-commutation relations and the orthonormal set of basis 
functions spans the discretized Hilbert space S^d- (Here, no charge conjugation or 
renormalization is needed as we will physically deal with finite systems only.) 



3.1. Spectrum 



The free part Hq of this Hamiltonian, i.e., without the external potential = 0, can 
be explicitly diagonalized. Performing a discrete Fourier transform on the lattice 



d{p) := J2 e-'^'-'^a^n, b{p) := ^ e-^^'^+'^H,^^,, 

n n 

for p G [—it/2£, +7t/2£), where the anti-commutation relations (7) imply 



(14) 



{a(p), a(g)t} = ^6{p - q), {a(p)t, a(g)t} = {a(p), a(g)} = 0, (15) 



we obtain 



Hn = - 



TT 



dp 



M fa(p)^a(p) - b{pyb{p)) + - cos{ip) fa(p)'^fe(p) + h{p)^d{p) 



.(16) 



This Hamiltonian can be diagonalized via a unitary transformation mixing the two types 
of operators 



A{p) 

Hp) 

with the explicit form 



Uip) 



d{p) 

Hp) 



(17) 



Uip) 



what leads to 



VE + M yjE-M 
^ \ -VE - M VE + M 



IT 



Tv/2i 



dpE{p) A{pyA{p) - B{p)^B{p) 



with the energy-momentum relation 



E{p) = WM2 + -cos2(£p). 



(19) 



(20) 
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Due to two effective types of fermionic excitations, A{p) and B{p), which enter 
the Hamiltonian with opposite energy signs, we obtain two symmetric energy bands 
separated by a gap of 2M. Each approximates the relativistic energy-momentum relation 
at the edge of the Brillouin zone, for p ^ ±7r/(2£). In order to obtain a positive 
Hamiltonian, we can perform the usual re-definition B{py <H- B{p) which corresponds 
to changing the vacuum state by filling all B[p) states with fermions. This is analogous 
to the Dirac sea picture in full quantum electrodynamics. In terms of this analogy, 
or A{p) create or annihilate an electron whereas B{p) or B{py create or annihilate a 
positron. An additional potential $„, if sufficiently localized in space, will not modify 
this spectrum but may introduce bound states (isolated eigenvalues) with energies lying 
in the gap [30]. 



3.2. Supercritical potential 



As an example for demonstration of supercriticality in the discretized system we consider 
the attractive Woods-Saxon potential^ 

^(^) = -TT5ra' ^'«'^>0' (21) 

for which the one-dimensional Dirac equation is analytically solvable in terms of 
hypergeometric functions. The bound-state energies E^. are determined by the equation 

B{-2g,g + s-XY _ js-gf-X^ _ r(x)r(y) 

Bi2g,-g + s + Xy {s + gy-X^' T{x + y) 

withs := y/M"^ — E'^/a, g := i-y/ {E + Wy — M'^/a, X := iW/a and depend continuously 
on the parameters of the potential [31]. 





-0.5- 



2 .^r^.- 



EA 



jpper band 



bound state 



2M 



lower band 



Figure 3. Left: The continuous Woods-Saxon potential ^(x) (red daslied line) and 
tlie discretized effective potential $„ -I- (— 1)"M (blue dots, connected with dotted line 
for visualization only). Right: The corresponding energy spectrum. (Not to scale.) 



Below, we compare the spectra of the continuous and the discretized Dirac equations 
with the Woods-Saxon potential. In the latter case, the discretized Woods-Saxon 
potential is defined by $„ := $(a;„) with Xn = i ■ n (see Fig. 3). For both cases 
we calculate numerically the lowest-lying bound state Eq as a function of the parameter 



I However, the discussed phenomena are generic and do not depend on the details of the potential. 



Optical lattice quantum simulator for QED in strong external fields 



9 



W (a and L fixed) which is a monotone function with dE^/dW < as long as 
— M < Eq < M . The dependence of the bound state energies Eq{W) for the Hubbard 
Hamiltonian (10) on the parameter W is quahtatively the same and quantitatively in 
very good agreement with the curve Eq{W) obtained as solution of (22) in the continuous 
case. At almost the same critical value Wcr ~ 2.878 both bound states disappear from 
the spectrum and turn into complex resonances§. The curves are compared in Fig. 4. 




Figure 4. Comparison of the bound state energy E as a. function of the potential 
strength W for the continuous (black dots and line) and discretized (blue dots and line) 
Dirac equation with Woods-Saxon potential (21) (a = 10, L = 1). On the right, the 
parabolic approach, E{W) ~ ±M -\-C±{W — W±Y , to the upper and lower continuum. 




Figure 5. Typical profile of a bound state (here E — — 0.37A/). Plotted are two 
components of the wave function '^]-^, 4'^ (blue dots and red diamonds) against x — n-d, 
supplemented by exponential asymptotes (blue and red solid lines). One function has 
always one zero and the other has no zeros like in the continuous Dirac equation. 

§ In one dimension this process is slightly more complicated than the well-known diving of the bound 
state into the continuum in 3 dimensions. Here, the bound state first turns at the threshold E ~ —M 
into an anti-bound state {E > — M and V AP — E'^ changes sign) and moves slightly up to eventually 
turn again and dive in the negative continuum E < —M as a resonance, cf. [32]. 
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4. Quantum simulator 

4-1. Ultra-cold atoms in optical lattices 

The main goal of this work is to propose a physical quantum system, composed of 
ultra-cold atoms moving in a specially designed periodic potential, which can effectively 
be described with a Fermi-Hubbard-type Hamiltonian of the form (10) and pseudo- 
relativistic dispersion relation (20) approximating the relativistic formula E"^ = M'^ +p'^ 
for energies around the Fermi-level. Moreover, excitations of the ground state should 
behave like particles and antiparticles and obey the Fermi statistics. 

Ultra-cold atoms loaded into optical lattices are conveniently described by effective 
discrete Hubbard-type Hamiltonians [33]. That kind of approximation is based on a 
construction of a set of orthonormal wave functions tpn ( Wannier functions) localized 
around the local minima of the potential W{x), giving rise to a regular grid of sites, 
and on the assumption that the single-particle Hamiltonian "H is approximately tri- 
diagonal in that basis, i.e. {ipn\'H\iprn) ~ for \n — m\ > 1. In consequence, the 
many-body Hamiltonian can be written in the form (10) with J„ = {ipn\'H\ipn+i) and 
Vn = {'ipn\'H\tpn) ■ There is a deeper connection between that approximation, in which 
only the lowest energy band is taken into account in the construction of the Wannier 
functions, and a spatial discretization of the theory in which the discretization step 
(equal to the period of the potential) introduces a natural cut-off in energies. In the 
latter approximation, the coefficients J„ and Vn correspond to the discretized kinetic 
(Laplacian) and potential terms in the Hamiltonian. In both approaches the energy 
spectrum is reduced to a single band. 

4-2. Bi-chromatic optical lattice 

It turns out that the emergence of a pseudo- relativistic dispersion relation, as in Eq. (20), 
is a quite universal phenomenon, see also [34, 35]. Imagine, we start with a periodic 
potential in one spatial dimension and introduce a small perturbation which breaks the 
original periodicity and is only periodic with the double period. This implies that the 
Brillouin zone [— 7r/(2£), 7r/(2£)] shrinks by a factor of two and that the lowest band 
splits into two sub-bands. Since, at the same time, the perturbation is small the energy- 
momentum relation E{p) at any given momentum p can only change by a small amount. 
In consequence, the perturbation will induce significant changes only in the vicinity of 
the momenta Pq = ±7r/(2£), i.e. edges of the shrinked Brillouin zone, at which it 
generates a small gap in the spectrum which separates the two branches of E{p), see 
Figure 6. For small perturbations, this gap will be proportional to the amplitude of 
the perturbation [34, 35]. Altogether, we reproduce the pseudo-relativistic dispersion 
relation in the vicinity of those points po- 

The conditions for a quantum simulator formulated above can be satisfied, in a good 
approximation, with ultra-cold fermionic atoms loaded into a one-dimensional optical 
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Figure 6. Sketch of the typical dispersion relation for a periodic potential (solid blue 
line) and perturbed doubly-periodic one (dashed red curves). 

lattice with the doubly-periodic potential 

W{x) = Wo sm'^(2kx) + AW sm^{kx) , (23) 
where k = 7r/(2£), by taking Wq ^ AW (see Fig. 7). Potentials of that form can be 




Figure 7. Left: Doubly-periodic potential W{x). Right: The corresponding energy 
spectrum with the gap 2Af « AW between the lowest two bands. (Not to scale.) 

obtained by superposition of two lattice-generating standing laser waves with different 
frequencies. Similar settings have already been obtained experimentally [36]. 

Unfortunately, it is not possible to find a closed analytic formula for the energy- 
momentum dependence in that potential. However, by applying a version of the WKB 
method for periodic potentials [37] to the doubly-periodic case we were able to derive [35] 
a spectral condition from which an approximate dispersion relation can be calculated 
analytically. That condition reads 

cos2(<l>/2) = (1 - T) sin2(A$) + T cos^iip). (24) 

where T{E) is the WKB-transmission coefficient through a single potential barrier 
[around a maximum of W^(a;)], $ := $i $2 and A<l> = $1 - $2 with the WKB 
phases 

^i{E) := f ' y^E - W{x) dx (25) 
Jy^ 

calculated between two consecutive turning points yi, Zi corresponding to the same 
potential minimum for which W{yi) = W{zi) = E and W{x) < E for yi < x < Zi. 
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The index 2 = 1,2 refers to two different types of the potential minima (lower and 
upper). 

For large Wq ^ AW, the lowest energy band is narrow and lies well below the 
potential maximum Wq. It implies small tunneling probability T{E) and small A$(i?) 
which are both relatively insensitive to E. The average phase can be approximated by 
a linear function $(-E) ~ a{E — Eq) around the value tt [first minimum of cos^($/2)] 
which leads to the effective universal relation 



E{p) -Eo^ ±^/WTJ^co^W) (26) 

where M := (1 — T) sin^(A$)/Q;^ and := T /a^. The approximation holds uniformly 
for all p G (— 7r/2£, 7r/2£) as long as AW <^ Wq (for more details, see [35]). Using again 
the WKB approximation, we can estimate the parameters 



^ ' \ AVE 



J.lV^exp , (27) 



where Er = k"^/ (2Matom) = tt^/ (SMatomO^) is the recoil energy and Matom the mass of 
the atoms moving in the doubly-periodic potential W{x). 

4-3. Wannier functions and sites 

Let us discuss the transition from the simply periodic potential to the doubly periodic 
one on the level of the associated Hamiltonian. Starting with the single-particle 
Schrodinger Hamiltonian describing atoms in an external potential 

H = -TTT^V^ + Wix) (28) 

and performing the standard steps we obtain, for the original periodic potential 
W{x) = W^o sin^(2A;x), the usual second-quantized Hamiltonian in momentum space 

-f^original = 7^ / dp Eonginax{,P) i'\p)^{p) , (29) 

27r 

where we consider the lowest band only with -E'originai(p) ~ —Jcos{ip) (for convenience 
we shifted the energy scale by the constant Eq what has no physical consequences). 
After switching on the doubly-periodic perturbation AW sin^^kx), the energy spectrum 
undergoes a transition -Eoriginai(p) E{p) ^ i^/M^ + cos^(£p) and the Hamiltonian 
becomes 



7r/2£ 



t 



with x{p) := ilj{p + tt) [assuming periodicity ipip + 2ti) = ip^p)] and 

\ + J2 COs2(£p) J- ^ ' 

This Hamiltonian has the same form as the one for the discretized Dirac equation (19) 
when we set J = 
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For the two separated energy bands there exist two separate sets of Wannier 
functions on the lattice: the "lower" and the "upper" centered at even and odd sites, 
respectively. But these Wannier functions turn out to be poorly localized on the lattice 
(somewhat analogously to continuous quantum field theory where free particles with 
fixed energy are not localized in space). In order to achieve optimal localization it is 
preferred to switch to the set of operators introduced already in (17)-(18) 

what transforms the matrix Ad via the similarity transformation Ai' = U'^M.U to 

^ - y jcosiip) -M ) ■ ^^^^ 

Now, going from the momentum to the site representation by inverting the Fourier 
transformation (14) 

dp e^'-'^dip) , = - / dp e^(2"+i)^^6(p) , (34) 



a2r, 

^ J -IT 121. ^ J-tt/2(. 

we obtain the free Hubbard Hamiltonian (9) with J = l/{2i). By this construction 
dl^ and &2n+i create two types of particles in two types of Wannier states exponentially 
localized at even and odd sites, respectively. However, they do not give rise to "positive" 
and "negative energy sites" as they mix energies from both bands. This can be best 
seen in the limiting case M <^ J where the Wannier functions [up to terms O [Mj J)] 

a2n = "^(^2n " ^2n), ^2n+\ = -^iX2n+l + ^2n+l) (35) 

are build from the difference and sum of the single-band Wannier functions for the lower 
and upper bands defined as 

■>TT/2e p /.7r/2£ 



p pTT/-zi p pixril 

- dpe'-'^ijip), Xn:=- dpe^-'^xip)- (36) 

^ J-TT/2e ^ J-TT/2e 



4 ■4- Physical parameters 

In order to discuss the experimental realizability of our quantum simulator, let us 
summarize the conditions on the involved parameters. Strictly speaking, the WKB 
approximation used above requires Wq 3> which then implies Er ^ J via Eq. (27). 
However, even if we relax these conditions to 

Wo>Er>J, (37) 

we still get qualitatively the same picture. What is crucial, however, is the applicability 
of the single-band Fermi- Hubbard Hamiltonian (10). To ensure this, we demand that 
the local oscillator frequency Wosc in the potential minima be much smaller than J. In 
addition, the continuum limit - i.e., that the discretized expression (6) provides a good 
approximation - requires J 3> M, i.e., l/£ 3> M. For the same reason, the change 
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A$„ = $„+i — $„ of the analogue of the electrostatic potential from one site to the 
next should be smaller than M. Over many sites, however, this change can well exceed 
the mass gap 2M, which is basically one of the conditions for the Sauter-Schwinger 
effect to occur. Finally, the effective temperature T should be well below the mass gap 
2M in order to avoid thermal excitations. In summary, the analogue of the e"'"e~ pair 
creation can be simulated if the involved scales obey the hierarchy 

Wo.c>J>M>T. (38) 

Let us insert some example parameters. The recoil energy Epi of ^Li atoms in an optical 
lattice made of light with a wavelength of 500 nm is around Eji ^ 7 fiK. Thus, if we 
adjust the potential strength to be Wq = 10 fiK, the hopping rate J would be around 
5 fiK which is still sufficiently below the local oscillator frequency Uqsc of around 34 fiK. 
Then a perturbation of AW = 1 fiK created by light with a wavelength of 1000 nm 
would induce an effective mass M of 500 nK and thus the effective temperature should 
be below that value - which is not beyond present experimental capabilities. 

5. Spontaneous pair creation on the lattice 

The above established analogy between the (discretized) second quantized Dirac field 
describing electrons and positrons in an electric field, on the one hand, and the (Fermi 
or Bose, see Sec. 5.2) Hubbard model describing ultra-cold atoms in an optical lattice, 
on the other hand, enables laboratory simulations of some of the relativistic phenomena 
of strong-field QED. The original Sauter-Schwinger effect [5] with a constant electric 
field E would correspond to a static tilted optical lattice with $(x) = Ex (the so called 
Wannier-Stark ladder, see e.g. [38]). For nonzero values of E > one would expect a 
constant rate of spontaneously generated particles and holes (cf. formula (1)), depending 
non-perturbatively on E. However, a constant electric field E is unrealistic from an 
experimental point of view. An electric field which is localized in space and time is 
simpler to handle both, experimentally and conceptually (see e.g. [6, 39]). Therefore, in 
the following we discuss the process of analogue spontaneous pair creation in presence of 
an external localized potential which will be slowly switched on to a supercritical value - 
admitting one bound state to dive into the negative continuum - and then switched off, 
as discussed in Sec 2. In presence of the attractive potential a bound state will form in 
the gap 2M between the two lowest bands (formed from the lowest band splitted due to 
the doubly-periodic perturbation). During the time-dependent process, the bound state 
will slowly reach the lower band and then turn into a resonance lying within this band 
(see Fig. 8). The resonance will then decay causing an instability of the Fermi state 
(our analogue vacuum state) which will spontaneously decay to an energetically more 
favorable state with a particle- hole pair present. The "particle" (an atom excited above 
the Fermi level) will stay bound by the attractive potential while the "antiparticle" (hole 
in the Fermi sea) will be in a scattering state. 
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Figure 8. Time-dependent bound state/resonance interpolating between two halves 
of the lowest energy band after splitting it by a doubly-periodic perturbation of the 
optical potential. (Bands and gap not to scale.) 



5.1. Experimental procedure 

Simulation of the "pair production" in the optical lattice requires preparation of the 
initial quantum state of the atoms corresponding to the Dirac sea in QED. This can 
be achieved by keeping a large value of /SW (separation of two bands by a large gap) 
during the cooling phase and generating a clear Fermi level at half filling of the lattice 
(all particles in the lower sites, see Fig. 9, top). Next, the value of AIV should be 
adiabatically decreased to values well below Wq to achieve M <^ J and allow the atoms 
for dispersion across the lattice. The atoms become delocalized but still the lower band 
is fully filled while the upper band remains empty (second picture in Fig. 9). Then 
the "external" potential $ (which mimics the electric potential) can be slowly switched 
on to reach a supercritical value. In that phase, the ground state ("vacuum") will get 
rearranged via tunneling from the lower band to the upper band (analogue of the Sauter- 
Schwinger effect, third picture in Fig. 9). Such an created particle-hole pair will tend 
to separate on the lattice so that when the potential is slowly switched off after some 
delay the pair will not be able to annihilate any more (fourth picture in Fig. 9). That 
mimics the well known spontaneous pair creation known from QED. Finally, in order to 
detect the "pair" in experiment, the value of AW can be adiabatically increased again 
thus leading to energetic separation of the created particle and hole represented by an 
atom in one of the upper minima and a missing atom in one of the lower minima (fifth 
picture in Fig. 9). The atom in one of the upper sites and the hole (missing atom) in one 
of the lower minima could be detected via site- resolved imaging [40]. Another option 
could be blue-sideband-detuned optical transitions which are resonant to the oscillation 
frequency Uqsc of the upper minima but not to those in the lower sites. 
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Figure 9. Sketch (not to scale) of the stages of the simulation (from top to bottom). 
Left: The solid red curves represent the sum of the optical lattice potential W{x) and 
the dynamically switched supercritical potential $(a;) as a function of position x while 
the dashed blue curves correspond to the effective electric potential The blue solid 
dots are particles and the empty circle is a hole. Right: Band structure and a bound 
state changing in time. 

5.2. Bose-Fermi mapping 

Since it is typically easier to cool down bosonic than fermionic atoms, let us discuss an 
alternative realization based on bosons in an optical lattice. To this end, we start with 
the Bose-Hubbard Hamiltonian 

H = -jJ2 [dUidn + dldn+i] + J2 ^-^ndn + ^ $^(4)' dl , (39) 

n n n 

which has the same form as the Fermi-Hubbard Hamiltonian (9) after replacing the 
fermionic c„ by bosonic dn operators, but with an additional on-site repulsion term U . 
For large U ^ J (which can be controlled by an external magnetic field via a Feshbach 
resonance, for example), we obtain the bosonic analogue of "Pauli blocking", i.e., at 
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most one particle can occupy each site dnl"^) ~ 0. Neglecting all states with double 
or higher occupancy, we can map these bosons exactly onto fermions in one spatial 
dimension via 




Via this transformation, the bosonic commutation relations [d^, d]^] = 5nm and 
[d\, d]^ = [dn, dm] = are exactly mapped onto the fermionic anti-commutation 
relations {c„, cj„} = 6nm and Cm} = {c|^, cj„} = 0. As a result, we obtain the 
same physics as described by the Fermi-Hubbard Hamiltonian (9). 



5.3. Interactions 

Apart from investigating the pair creation probability for space-time dependent electric 
fields E{t, x), this quantum simulator for the Sauter-Schwinger effect could provide some 
insight into the impact of interactions. For example, including dipolar interactions of 
the atoms, we would get the coupling Hamiltonian DnmCncl^CnCm with Dnm oc \n — m\~^. 
As an example for permanent dipole moments, we may consider ^^Cr atoms possessing a 
rather large magnetic moment. However, the associated interaction energy Dnm would 
be below one nano-Kelvin and thus probably too small to generate significant effects. 
Therefore, let us consider externally induced dipole moments. For example, ^Li atoms 
can be electrically polarized by an external electric field of order 10® V/m (which can 
be realized experimentally) such that the induced electric dipole moments generate 
interaction energies up to a few /iK. By aligning the atomic dipole moments parallel 
or perpendicular to the lattice, we may even switch between attractive Dnm < and 
repulsive Dnm > interactions. Note that this goes far beyond the simulation of the 
classical Dirac equation and requires the full quantum many-particle Hamiltonian. 



Appendix A. Potential localized at one site (delta-like) 



In order to compare the discretized Dirac Hamiltonian with its continuum version, we 
consider an example where both can be solved analytically. This is possible for a Dirac 
delta like potential localized at one lattice site 

1 1 ^0,n 



with 



(p{p) = (f) 



(A.l) 



which corresponds to 0(x) = (j)6{x) in the continuous case £ — )■ 0. It has the advantage 
that a closed formula for the bound-state energy can be found analytically. The Hubbard 
Hamiltonian (16) with added potential reads then 



H = - 

vr 



+ 



■K/2e 



dp M (a(p)^a(p) — b{pyb{p)] + J cos{ip) id{pyb{p) + b{pyd{p) 



7r/2£ 



dpi a) (pi 



7r/2£ 



TT/2e 



dp2 d{p2) + 



■k/21 



■k/21 



dpi &"^(pi) 



Tr/2e 



dp2 b{p2 



+ 



(A-2) 



-7r/2£ 
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where we set J = Obviously, it satisfies H\Q) = for the vacuum vector |0) such 
that a(p)|0) = h{p)\Q) = 0. However, we want to find a one-particle eigenstate |x) of the 
Hamiltonian H with eigenvalue A lying in the gap —M < X < M (which corresponds to a 
bound state for the discretized Dirac equation). Then we want to study the dependence 
of the eigenvalue A on the strength of the potential 0. The general one-particle state 
can be written as 

dpA{p)a\p) H — / dp B[p)h\p)dp 



\x) 



TT 



|0) (A.3) 



'-it/21 " J-7r/2£ 

with two complex functions A[p) and B{p). By projecting the eigenvalue equation 
H\x) = \\x) oil the states (0|a(p) and (0|6(p) we arrive at a system of equations 

M-A, Jcos(£p) A / A{p) \ _ f^'^^ 
Jcos{£p), -M-A j I B{p) 




(A.4) 



= :M 

Integration of both sides over p leads to the relations 

7r/2£ 



^ ^j:if^^dpJcosiep)B{p) 



A = - A{p)dp 



'-./2E A-M 



r-/2^ ^ r^'^,i,dpJcos(ip)A(p) 

B:=^l B{p)dp = -^-''[ \^ . (A.5) 



J-^/2i A - M - ^ 

Observe that the right-hand side of (A.4) is equal to —(j){A^ B)^ . We can invert the 
matrix A^, whose determinant det = A^ — — Jcos^{ip) < never vanishes, to 
obtain 




.(A.6) 



A{p) \ / A + M, Jcos{ip) 

B{p) J ~ A2 - M2 - Jcos2(£p) 1^ Jcos(£p), A-M 

Integrating both sides over p again and using 

i_ r'^^ dp 1 

^ J-^/2l A2-M2- JcOs2(£p) ~ ~VM2- AVM2 + J2- A2 ' 

^/^^ j^osmdp ^^^^^ 



(A.8) 



J-n/2e A2 -M2 - Jcos2(£p) 
we obtain consistency conditions 




(A + M)A 

VM2 - A2VM2 + J2 _ A2 I (A - M)B 



which cannot be satisfied at the same time except when (at least) one of the constants 
A, B vanishes. Assume first, it is 5 = (the case A = will be discussed below). Then 
we need to solve the algebraic equation 

A + M 



VM2 - A2VM2 + J2 - A2 



(A.9) 
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Figure Al. Left: The eigenvalue X/M as a function of the potential's strength (f) 
for the discretized Dirac equation (red solid line) compared to the continuous Dirac 
equation with the delta potential (blue dashed line). 

Right: The momentum distribution A{p),B(j)) for a bound state with A = 0.8A/. 



for the function A(0). It has, in general, three solutions what can be better seen from 
(M-A)(M=^ + J^-A^) = (/)2(M + A). (A.IO) 



For = 0, these three solutions start from the points M, i-y/M^ + J^, i.e., the edges 
of the continuous spectrum. We are interested in the perturbations of the eigenvalue 
A = M for negative values of 0, i.e. for a bound state separating from the bottom of the 
upper band. For small 0, A(0) behaves like A(0) = M(l — 20^). For increasing values 
of 1 01 it monotonically decreases to — M but never reaches this value (more precisely, 
A(0) = — M+2M/0^ for |0| 3> 1) - see Fig. Al. It shows that there is no supercriticality 
in this potential, i.e., crossing of the value A(0) = — M at finite 0, what is analogous to 
the continuous case 0(a;) = (j)6{x) in which A(0) = M(l — 0^)/(l + 0^) [41]. 
From equation (A. 6) we can also find the eigenvector to the eigenvalue A 




A X + M 



A2-M2- Jcos2(£p) \ Jcos{ep) 



(A.ll) 



where the value of A is to be determined from the normalization condition 

"271 



1 [ \p {\A{p)\' + \B{p)\') = 1 . (A.12) 



27r 

Since both |y4(p)p and \B{p)\'^ are even functions of p, we have 

^I'^dp {\A{p)\'+\B{p)f) Jcosiip) = 0, (A.13) 

which refiects the fact that the (discretized) momentum id/dx — )■ J cos{ip) vanishes in 
the bound state. Note also that the condition B = implies Bq = what means that 
the "antiparticles" are repelled from the site at which the potential is localized. 

The case A = is fully analogous and the solutions can be obtained by a symmetry 
transformation: A — t- —A and — > —0. The bound state emerges then from the lower 
band (negative "continuum" ) at A = —M and goes up for positive potentials 0. 
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Remark 1. The considered potential is localized at one site n = in that sense that 
it interacts with both types of particles via and Bq. But in fact, and Bq are two 
different (neighbouring) sites brought to = by a convenient renumbering. It is also 
possible to consider the potential to be localized in such a way that it interacts with 
only one type of particles, say via ao- Then the equations get slightly modified (the term 
13 disappears at some places) and we obtain 5 = as consequence of (A. 6) and (A. 7). 
From that point on, the solution is identically the same to the previous case. It means 
that it plays no role whether we consider potentials localized at one site interacting with 
only one or with both types of particles because, in the latter case, the solutions split 
into two symmetric cases of the former type. 

Remark 2. The above method of calculating A((/)) works only for potentials for which 
^ Jq"" dp (pip— q) det A^(p)~^ is independent of g, i.e. for 0(p) = const which corresponds 
to (h n ~ ^n,o- Unfortunately, it cannot be generalized to more complex potentials in a 
simple way. 
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